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We study the magnetic properties of metallic Gai-j^Mn^^As. We calculate the effective RKKY 
interaction between Mn spins using several realistic models for the valence band structure of GaAs. 
We also study the effect of positional disorder of the Mn on the magnetic properties. We find 
that the interaction between two Mn spins is anisotropic due to spin-orbit coupling both within 
the so-called spherical approximation and in the more realistic six band model. The spherical 
approximation strongly overestimates this anisotropy, especially for short distances between Mn 
ions. Using the obtained effective Hamiltonian we carry out Monte Carlo simulations of finite 
and zero temperature magnetization and find that, due to orientational frustration of the spins, 
non-collinear states appear in both valence band approximations for disordered, uncorrelated Mn 
impurities in the small concentration regime. Introducing correlations among the substitutional Mn 
positions or increasing the Mn concentration leads to an increase in the remnant magnetization at 
zero temperature and an almost fully polarized ferromagnetic state. 

PACS numbers: 75.30-m,75.30.Cr,75.30.Hx,75.50.Pp 
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I. INTRODUCTION 

Recently, there has been a surge of interest in 
the more than 30 year old field of diluted magnetic 
semiconductorsi that has been largely motivated by 
potential application of these materials in spin-based 
computationSiMi^ devices. In particular, the discovery of 
ferromagnetism in low-temperature molecular beam epi- 
taxy (MBE) grown Gai-jjMuj^As has generated renewed 
interest. In this material Curie temperatures as high as 
Tc « 160K have been observedi^ 

It is typically rather difficult to dissolve magnetic 
impurities in a semiconductor using conventional MBE 
techniques, because the magnetic ions usually tend to 
segregate from the host. As a result, these convention- 
ally produced magnetic semiconductors often exhibit spin 
glass physics at low temperatures and other, undesir- 
able, "clustering" effects. Recently Ohno achieved the 
key breakthrough^ by growing Gai_2:Mn2:As at close to 
room temperature, which minimized Mn positional re- 
laxation during the growth process, and resulted in a 
material with intrinsic ferromagnetic properties, with a 
Curie temperature Tc ^ llOK. 

Post-growth annealing of the samples was shown 
to induce changes both in the lattice and positional 
defectsSiSiSiiS as well as in the hole concentration i!^ 
Properties such as magnetization and Curie tempera- 
ture, were found to depend very sensitively on details 
of the post-growth annealing protocol2*i2ii^*i^. In par- 
ticular, the measured T — temperature magnetization 
has been found to be much less than one expected based 
on the nominal concentration of Mn ions.^- The remnant 
magnetization and the Curie temperature Tq were both 
found to increase upon annealing, but the magnetization 



reached only about half of the expected value even in 
the case of optimally annealed samplesSiiLii. Accord- 
ing to recent experimental results, much of the missing 
magnetization is probably due to interstitial Mn ions, 
which compensate valence band holesjSiS and, presum- 
ably, also bind antiferromagnetically to substitutional 
Mn ionsii^ As a result, the nominal concentration of Mn 
ions, X, is usually substantially larger than the concen- 
tration of active Mn ions (participating in the ferromag- 
netism), Xactivo, which is larger than the concentration 
of mobile holes in the valence band (or possibly impurity 
bandiSiSfl), c = a:active /, with / < 1 the hole fraction. 
[In the remainder of this paper - unless explicitly noted 
otherwise - we present the results in terms of the active 
Mn concentration: x = Xactivo-] However, especially for 
samples with a lower doping level and/or shorter anneal- 
ing times, even this substantial difference in active and 
nominal Mn concentration seems to be insufficient to ex- 
plain the total amount of lost of magnetization: In some 
of these samples the remnant magnetization can be sub- 
stantially increased by a relatively small magnetic field, 
clearly hinting at a non-collinear magnetic statciii^ 

Many of the properties of Gai-a^Mn^^As depend cru- 
cially on defects, and understanding these defects is 
currently a subject of intense research 
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Gai_2,Mn2,As is a complicated material: Besides contain- 
ing many defects of unknown origin, it is very disordered, 
close to a localization transition with a mean free path 
of the order of the Mn-Mn distance, has a complicated 
band structure, has strong spin-orbit coupling and has a 
large exchange coupling between the localized Mn spins 
and the valence holes. It is therefore almost impossible 
to study this material without making some approxima- 
tions. 



As pointed out in Ref. [23, spin-orbit coupling can 
induce frustration and non-collinear states in a disor- 
dered magnet. In this paper, we study the effect of 
spin-orbit couphng in the metaUic Hmit (i.e., holes are 
assumed to reside entirely in the valence bands of GaAs) 
of Gai-xMn^As, and study how it may influence the 
structure of the ferromagnetic state. To this purpose, we 
will follow the route outlined in Ref. [23 and compute 
the effective interaction between the spin of two substi- 
tutional Mn impurities by doing perturbation theory in 
the exchange interaction (RKKY interaction). It is im- 
portant to note that other mechanisms also exist that 
could lead to a non-collinear statei^SiSi Recently, several 
studies of RKKY models of Mn spin-spin interactions 
in GaMnA a^^'^i^?i^?i?9'?-^i?? have been reported, but these 
studies have not investigated in detail the effects of posi- 
tional disorder on the magnetic properties using realistic 
band structure models. 

In this work, we neglect the effect of other defects 
and the strong scattering potential on the charged sub- 
stitutional Mn impurities. This approximation - which 
is quite common in the literature^SiS - is expected to 
provide qualitatively reliable results relatively deep in 
the metallic state, far away from the localization tran- 
sition, where screening is expected to reduce the effets 
of charged impurities. Quantitative agreement with ex- 
periment will require treating both charging effects and 
defect correlations in a realistic way^i*^ Furthermore, 
scattering from charged impurities becomes even more 
important in the dilute limit, and is ultimately respon- 
sible for the localization transition that occurs at small 
concentrations. In this dilute limit, the strong potential 
scattering off Mn impurities can be treated nonpertur- 
batively, and results similar to those reported here are 
obtained.i2i2^ 

In the absence of disorder the top of the valence band of 
GaAs can be described within the framework oi k-p per- 
turbation theory, which also accounts for the spin-orbit 
coupling in this material, and gives a good description 
of the band structure around the F point. In particu- 
lar, we study two different forms of the k -p perturbation 
theory. First we study the effective spin-spin interaction 
within a simplified version of the four band model, the 
so-called spherical model, ^® where only the topmost four 
spin j = 3/2 bands are kept and the distortion of the 
spherical Fermi surface is neglected. In this case ana- 
lytical results can be obtained for the effective Mn spin- 
spin interaction, ^5 Then we study the Mn spin-spin in- 
teraction using a six band modelSiSi (Kohn-Luttinger 
IIamiltoniari2i) , where the spin-orbit split j = 1/2 band 
is also taken into account. In the latter case it is not pos- 
sible to evaluate the interaction kernel analytically, and 
one must resort to numerics. 

Although we find in both calculations a strong spin- 
orbit induced anisotropy in the spin-spin interaction for 
typical Mn-Mn distances, the spherical model largely 
overestimates the size of the anisotropy for small dis- 
tances. We find it especially instructive to discuss this 



result in comparison with those published recently in a 
nice and intriguing paper by Brey and Gomez-Santos^ 
While our result agrees qualitatively with the one ob- 
tained in Ref. 39] , on a quantitative level it is completely 
different. Brey and Gomez-Santos estimated that the 
largest value of the anisotropy is of the order of ^ 10^^. 
In contrast, in the present calculation we find that the 
smallest value of the anisotropy is around 1% for nearest 
neighbor Mn ions at a distance of ~ 2.5 A, and it increases 
continuously with distance, until it reaches a value of the 
order of ~ 20 % for typical Mn-Mn distances. These 
numbers roughly agree with the results obtained in the 
very dilute limit for a Mn2 "molecule" within the six 
band model»^ 

There are basically two reasons for the discrepancy be- 
tween the two results: 

(1) Presumably to avoid convergence problems, Brey and 
Gomez-Santos introduced a short distance cut-off oq, and 
replaced the exchange interaction between the Mn core 
spins and valence holes by a non-local interaction. The 
cut-off scheme they introduced, however, is not com- 
patible with the general structure of the exchange in- 
teraction derived in the microscopic theory of quantum 
impurities.^*' As we discuss in detail in Sec. IIIII the cor- 
responding momentum cut-off introduced should not be 
smaller than ~ 2fci?, the typical momentum transfer dur- 
ing backscattering. Unfortunately, the cut-off oq = iA 
used in Ref. 39] did not satisfy this criterion, and as a 
consequence, the results of Ref. [33 showed a very strong 
dependence on ag: In particular, for somewhat smaller 
values of oq — 2.42A Brey and Gomez-Santos found an 
anisotropy of a few percent, which roughly agrees with 
the one we find (« 1 %) within a cut-off scheme com- 
patible with the general theory of magnetic interactions 
for nearest neighbor Mn ions. The strong dependence of 
the Brey- Gomez- Santos result on the cut-off parameter 
was also pointed out and questioned recently by Timm 
and MacDonald in Ref. '3l'| . We emphasize that with our 
cut-off scheme the anisotropy depends only weakly on the 
value of qq . 

(2) Secondly, in Ref. [33 it has been assumed that the 
anisotropy is largest for the shortest Mn-Mn separa- 
tions. However, for hole concentrations p < 3 nm~^ 
one actually expects on general grounds that the asymp- 
totic form of the RKKY interaction is reasonably well- 
approximated within the four band model, which as- 
sumes infinite spin-orbit splitting Aso^ a-nd that the 
anisotropy increases with Mn-Mn separation. This is in- 
deed supported by our numerical results: For the very 
small separations for which Brey and G'omez-Santos car- 
ried out their calculations we indeed find only a ~ 1% 
anisotropy, while for typical Mn-Mn separations we find 
a rather large anisotropy of the order of ^ 20%. 

Having the interaction kernels at hand, we considered 
a distribution of Mn ions and carried out finite temper- 
ature Monte Carlo simulations to characterize the mag- 
netic properties of the system. At present, it is unclear 
from the experiments to what extent the positions of sub- 



stitutional, magnetically active Mn are correlated during 
the growth process. In order to gain insight into how Mn 
positional correlations may affect the magnetic properties 
of a material like GaMnAs, we studied the dependence of 
the Curie temperature, the saturation magnetization, the 
shape of the magnetization curve and the ground state 
spin distribution on Mn positional correlations. These 
correlations were introduced by allowing repulsive inter- 
actions between the Mn ions and relaxing them through 
a zero temperature Monte Carlo Metropolis algorithm in 
a way similar to Ref. |3J|. We started the simulations 
from a set of completely uncorrelated initial Mn posi- 
tions, and allowed the substitutional Mn ions to move. 
We then studied the magnetic properties of the spin sys- 
tem as a function of the Mn positional relaxation time, 
including the completely uncorrelated configuration. 

Our main results are as follows: Within the spheri- 
cal approximation we find a strongly disordered (non- 
collinear) ferromagnetic ground state for unrelaxed Mn 
positions, where the ground state magnetization can be 
reduced by as much as ~ 50% with respect to its satura- 
tion value. Upon relaxing the Mn positions, however, the 
remnant magnetization increases, and for a fully relaxed 
lattice one may recover as much as 95% of the saturation 
value. In contrast to our earlier expectations,^^ the ob- 
tained non-coUinear state is not a spin-glass: Although 
a finite system displays hysteresis, the coercive field we 
compute decreases with system size, implying that the 
hysteresis vanishes in the infinite system size limit. The 
frustration effects found within the spherical model are 
very sensitive to the hole fraction / characterizing the 
degree of compensation, and become more pronounced 
for larger hole fractions. 

Spin-orbit coupling induced anisotropy does not have 
such a dramatic effect if we use the kernel obtained within 
the six-band model. For a substitutional Mn concentra- 
tion of a; = 3%, hole fraction / = 0.4, and unrelaxed Mn 
spins we find that the magnetization is reduced by ~ 5%, 
relative to the fully polarized state, corresponding to a 
typical non-coUinear angle of the order of 6* ~ 16 degrees. 
This non-coUinearity is, however, almost completely sup- 
pressed if we let the Mn positions relax using the Monte 
Carlo algorithm mentioned above. By increasing the sub- 
stitutional Mn concentration to a; = 5% the ground state 
becomes almost fully collinear even without relaxing the 
Mn positions. The results do not depend too strongly on 
the level of compensation (hole fraction) . Our results also 
indicate that anisotropy effects may be more important 
for smaller Mn concentrationsfi^, 

We emphasize again that the concentration x above 
denotes the concentration of active Mn ions, i.e. the 
concentration of those Mn ions that participate in the 
formation of the ferromagnetic state. Due to compen- 
sation effects induced by interstitial Mn impurities, x is 
usually much less than the nominal concentration of Mn 
ions. In fact, samples with a nominal concentration of 
5% may easily have an active Mn concentration in the 
range x < 0.02 — 0.03. In this concentration range the 



impurity band approach of Refs. [23 and [13I may be 
more appropriate. 

In agreement with the results of Ref. [41| , we also find 
that in all cases Tc decreases as we relax the Mn impu- 
rities, irrespective of the specific form of effective inter- 
action used, and M{T) becomes more mean-field like at 
the same time. 

This paper is organized as follows. In Sec. |n] we 
present the RKKY calculation within the Baldereschi- 
Lipari spherical approximation for the effective interac- 
tion between two Mn spins. Several of the most technical 
details and lengthy expressions have been relegated to 
Appendices 1^ and IBl Sec. lllll is devoted to the computa- 
tion of the effective Mn spin-spin interaction within the 
full six band model. In Sec. lIVI we discuss the results of 
classical Monte Carlo calculations for the magnetization, 
the ground state spin distributions, the Curie tempera- 
ture and the dependence of these quantities on positional 
disorder and hole concentration using the effective spin- 
spin interactions computed in SectionsHlland UTTI Finally, 
in Sec. we summarize our main conclusions. 



II. EFFECTIVE MN SPIN-SPIN INTERACTION 
IN THE SPHERICAL APPROXIMATION 

In the present section, we study the effective RKKY 
interaction between two Mn impurities, by using the so- 
called spherical approximation to describe the valence 
band hole fluid. The top of the valence band of GaAs con- 
sists of six p-bands. Four of these bands become fourfold 
degenerate at the F point and have spin j = 3/2 charac- 
ter, while the other two are split off by A50 ~ 0.34 eV 
and have spin j = 1/2 character i^ 

The spherical approximation consists of setting the 
spin-orbit splitting Aso to infinity, keeping only the spin 
j = 3/2 bands, neglecting the warping of the Fermi sur- 
face due to the cubic symmetry of the crystal, and ap- 
proximating the dispersion of the holes by the following 
expression!^ 
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Here J^p = \[jaJ0 + jpja) ~ ■|<5a/3Tr{j„J>} is the 
quadrupolar angular momentum of the holes. The linear 
momentum tensor of the holes Pap is defined likewise. 
The couphng v — (673 + 472)/57i « 0.77 characterizes 
the strength of the effective spin-orbit coupling within 
the i = 3/2 band, and the 7i's are the so-called Luttinger 
parameters'^, characterizing the band structure of GaAs. 
Physically the parameter v describes how the spin of a 
i — 3/2 valence hole couples to its orbital motion. 

The advantage of the spherical approximation, Eq. |^, 
is that it makes it possible to obtain analytical results 
while it gives a realistic description of the top of the va- 
lence band. As we will see, this latter fact is somewhat 



misleading: While providing a qualitatively correct de- 
scription of the valence bands of Gai_a;Mna;As, it gives 
quantitatively incorrect results for the RKKY interaction 
between two local spins. 

In GaMnAs subsitutional Mn ions are believed to be in 
a Mn^"*" configurarationiSii^ii^ thereby behaving as a neg- 
atively charged scatterer and having a core spin S = 5/2. 
This core spin then couples to the spin of the valence hole 
fluid through an exchange interaction that, for a single 
Mn ion, takes the following form within the spherical ap- 
proximation: 



i/i„t(R) = GS-j(R), 



(2) 



where j(R) denotes the valence holes' spin density, at the 
location of the Mn ion, R, and G > is an antiferromag- 
netic coupling. 

To proceed let us first show that the eigenstates of 
Hgp are chiral, i.e., that the spin is quantized along the 
direction of propagation. First we note that Hsp for holes 
propagating along the z— direction is easily diagonalized: 
p — (0,0, pz) so that 
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Thus, for holes moving in the z— direction, 
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with the plus and minus signs standing for jz = ±1/2 
and jz = ±3/2, respectively. This defines the light hole 
mass, mi = ^^ r!\,A ~ 0.07m and the heavy hole mass 
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0.5™ where m, is the bare electron mass. 



For holes propagating in a general direction, the wave 
function can be constructed by spin 3/2 rotations so that 
// = j • k = ±3/2 for heavy holes and ±1/2 for light holes 
propagating along direction k. 

The kinetic energy Hgp greatly simplifies in the basis of 
these chiral states, and takes the following form in second 
quantized formalism: 
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(4) 



Here cL denotes the creation operator of a hole with 

wavevector k and spin projection /i along k. Diagonal- 
izing Hsp, however, produces a strongly momentum de- 
pendent carrier-impurity interaction: 
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(5) 

Here J^^,(k,k') ^ E,,/ Dl(M)jf^'D,,,,ik') where i^(k) 
is the spin 3/2 rotation matrix and V is the volume of 
the sample. The rotation matricies are defined here in the 



usual way as Dj^(k) = [e~*~''~e~*~'i~]j^, where (j) and 9 
are the usual azimuthal and polar angle of the direction 
k= (sin0cos(/), sin 9 sin cj), cos 9). 

Having diagonalized Hsp , we now consider two Mn ions 
at positions R and the origin, do perturbation theory in 
the exchange coupling G and compute the RKKY inter- 
action between two Mn spins.—- The first order correction 
to the single particle eigenstate |k, /i)o of Hsp is given by 
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where the sum is over all states except for |k, /i)o, and 
Ek^ = ^^/Sti^. With these single particle states in hand, 
we can then compute the shift of the ground state energy 
as: 
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(7) 



where we have dropped a piece of the energy indepen- 
dent of relative spin orientations. This equation can be 
rewritten as 



ifRKKY = -^^"(R)X"''(R)5'3(0) 



(8) 
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with an obvious definition of the kernel K°'l^ . Within the 
spherical approximation this result further simplifies to 



Hcff = -ifpara(i?)Sll • S'l - A'pe,p(i?)S^ • S 



(9) 



where R = |Ri — R2I and S; is the Mn spin at posi- 
tion Ri. Here S-'-/" refers to the component of the spin 
perpendicular/parallel to the axis joining them. 

The computation of the kernels K\\ and K± is not 
straightforward even in the case of the spherical approx- 
imation. Their explicit expressions and some details 
about how to obtain them are given in Appendix ^ 
These expressions together with Eq. © constitute the 
main result of this section of the paper. 

For more physical transparency, it is worth defining 
the following rescaled kernel. 
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where we have defined the dimensionless coupling as 
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with m,h is the heavy hole 



V2772 ^F 

mass and ep is the Fermi energy measured relative to the 
"top" of the valence band. The positional dependence 
of Cpaia/pcrp is shown in Fig. ^ where we also plotted 
the contribution of the heavy holes, giving the dominant 
contribution to the interaction kernel. The rather in- 
volved explicit expressions for the dimensionless kernel 



C{y) are given in Appendix ^ The asymptotic forms 
for y = kp^hR — ?■ of the kernels C are given by 
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where m; and irih are the light and heavy hole mass 
defined just below Eq. (jSJ. Specifically, in the limit of 
m^ ^ oo we simply obtain C^^^^y « l/2y, C^-^ « 1/y. 
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FIG. 1: Effective RKKY kernel in the 4-band spherical ap- 
proximation. The contribution from the heavy hole sector is 
also shown. At typical Mn-Mn separations the structure of 
the kernel is very sensitive to the value of the hole fraction /. 
For / = 0.10 one has y*^P = d^^'^kp^h = 1.6, while for / = 0.25 
we have y^^^ = d}'^^kF,h = 2.2. In both cases the perpendic- 
ular part of the kernel is ferromagnetic, but for / = 0.25 the 
parallel component is antiferromagnetic and is comparable in 
absolute value to the perpendicular component. 

Eq. ©, was first obtained in Ref. [23 using only the 
heavy hole sector of the valence band structure. There 
it was shown that Eq. ^ leads to a spin-disordered 
ground state whose signatures seem to be present in some 
experiments ,iS4i 

At typical Mn-Mn separations the kernel is very sen- 
sitive to the value of the hole fraction. Hole fractions in 
the range / = 0.10 — 0.25 corresponds to y'-^'^/^o.io = 

dtyp^/=0.10^^ g ^^^ ytyP^^„25 = d'y^ k^f-^' ^2.2. In 

both cases the perpendicular part of the kernel is largest 
and is ferromagnetic. However, for / = 0.10 the paral- 
lel component is ferromagnetic while for / = 0.25 the 
parallel component is antiferromagnetic and of compara- 
ble strength to the perpendicular component. The rela- 
tive size of iiTpcrp and ifpara means that spins prefer to 
ferromagnetically align perpendicular to the axis joining 
them for small enough hole concentrations. In general 
this leads to frustration and non-coUinear magnetism. 



In Sec. IIVI we study the finite temperature behavior 
of the Mn spins using classical Monte Carlo techniques. 
There we shall assume that the spins behave as classical 
objects, and we replace Eq. © by 



H,s ^ -Jn{R)n^l -f^ll 



j±{R)ni -n^ 



(11) 



where now J7i and $12 denote the two unit vectors spec- 
ifying the direction of the spins, and the exchange cou- 
plings are given by 



II/-L ^ ^J (-^para/pcrp 

with the exchange energy defined as 



(12) 



(13) 



The precise value of Ej depends on the specific value 
of the unknown coupling G and the hole concentration, 
p, with which it scales as Ej ^ e|, ^ p^/'^ within the 
spherical approximation. 



III. EFFECTIVE MN SPIN-SPIN 
INTERACTION IN THE SIX BAND MODEL 

In the previous section we studied the effective RKKY 
interaction between two Mn impurities within the spher- 
ical approximation. In the present section we shall study 
the effective Mn spin-spin interaction by using a more 
realistic description of the top of the valence band, the 
so-called six band model, where the Hamiltonian reads 



H 
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^Lutt + ^int 7 

Jpd(s(Ri) ■ Si 



-S(R2)-S2), 



(14) 
(15) 



with T^Lutt being the so-called Kohn-Luttinger Hamilto- 
nian, and s(R) the spin density of the holes at position 
R. Note that s here is not a spin 3/2 spin object. It is a 
6x6 dimensional matrix that represents the spins (1/2) 
of the holes in the three p channels that constitute the 
top of the valence band (for an explicit definition see e.g. 
Ref. H). 

The six-band Hamiltonian approximates rather well 
the band structure and the spin content of the hole states 
in the vicinity of the T point of the Brillouin zone. How- 
ever, it is impossible to evaluate the RKKY interaction 
analytically using -ffLutt, and numerical calculations are 
needed. 

To perform the numerical calculations, we followed 
the procedure of Ref. 39] : We generated a mesh in 
momentum-space, computed the eigenenergies of -ffLutt 
and the corresponding eigenfunctions, and ordered them 
according to their energy. We then computed the matrix 
elements of Hint in this basis by treating the spin of the 
Mn ions classically, and fixing their direction. Next we 
diagonalized i?Lutt + ^int to obtain the single particle 
hole states for this orientation of the spins, while only 
keeping states below the cut-off energy, -Ecutoff ■ Finally, 



we computed the ground state energy of the whole sys- 
tem as a sum of the energy of the occupied single particle 
states, and thus obtained the ground state energy of the 
system for a fixed number of holes as a function of the 
direction of the two spins. 

The effective couplings have been extracted from the 
ground state energy in the same way as in Ref . 6^ : First 
we placed the two Mn ions in a distance R along the z- 
direction and computed the ground state energies for a 
fixed number of holes in a case where the two spins were 
(1) parallel and pointing along the z-direction (E^^), (2) 
oppositely aligned along the z-direction (E^i), and (3) 
parallel along the x direction (E^^). In this arrange- 
ment the effective coupling and anisotropy can be defined 
as: 
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TT 



(16) 
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There are, however, a few technical details that should 
be discussed. First, the fc-points in Brillouin space must 
be generated in such a way that they respect the cubic 
symmetry of the crystal, otherwise the obtained effec- 
tive interaction breaks the cubic symmetry. Due to the 
cubic symmetry, however, the electronic levels typically 
become extremely degenerate for Jpd = and for a typi- 
cal mesh of k points this degeneracy can reach numbers 
as high as 96 in a mesh size of 1000. One therefore ex- 
periences large but apparently systematic fluctuations of 
the exchange energy computed as one fills the single par- 
ticle energy levels for Jpd ^ one by one (see the inset of 
Fig. 12)), and results converge very slowly with increasing 
mesh size. There are, however, special points where the 
number of holes is such that for Jpd = an integer num- 
ber of degenerate energy "shells" are occupied. At these 
points the calculations seem to converge somewhat faster, 
and an accuracy as high as 10 — 20% can be reached. This 
technique has been exploited in Ref. 39] , and we will also 
use it to compute the effective interaction kernel. 

Another very important issue is the cut-off scheme 
used: To facilitate convergence, it is important to intro- 
duce a cut-off for the exchange coupling in momentum 
space. Brey and Gomez-Santos therefore replaced Jpd in 
Eq. (|15() by a non-local interaction between the Mn spin 
and the valence hole, Jpd — > Jpd(r — R-), and introduced 
a corresponding cut-off for the exchange interaction in 
k-space of the form: 



Jpdi'^j k j — J J 



pdi 



-(k-k') = a^/2 



(18) 



where ag is the range of the non-local interaction. 

Unfortunately, there is a serious problem with the cut- 
off scheme (I18|l for large values of aQ-. Eq. H18|) gives a cut- 
off for the momentum transfer. On the other hand, it is 
well-known from the Kondo literature that the exchange 
interaction deduced from the more fundamental Ander- 
son impurity problem has the following structures^ 



Jpd(k,k')-liFk* 



1 



U + ed 



(19) 



where e^ < is the energy of the d-level, U is the Hub- 
bard interaction, and the Vk's denote the k-dependent 
hybridization with states within the Brillouin zone. The 
Vic's do have a momentum cut-off for momenta of the 
order of the size of the Brillouin zone, however, J(k, k') 
clearly factorizes, and there is obviously no cut-off for 
the momentum transfer across the Fermi surface,'^^ 

Since back scattering with large momentum transfer is 
responsible for much of the leading term in the asymp- 
totic RKKY interaction, it is essential to use a cut-off 
scheme which does not influence these large momentum 
transfer electron-hole excitations in the vicinity of the 
Fermi surface, and is consistent with Eq. p9|l . In fact, the 
value flQ = 4 A chosen by Brey and Gomez- Santos seems 
to be too large: They flnd that the relative strength of the 
anisotropy is very sensitive to the precise value of oq , and 
increases dramatically as they decrease it to ag — 2.5A, 
suggesting that the results obtained by Brey and Gomez- 
Santos with oq — AA are not reliable. 

To avoid this problem we therefore used the following 
cut-off scheme. 



Jpd (if, k j — Jpd e 



\li 



(20) 



which is consistent with Eq. p9|l . and does not suppress 
back-scattering. Using this latter cut-off scheme and 
flo = 4 A we flnd a magnetic anisotropy more than one 
order of magnitude larger than was reported in Ref. |39j | . 
Furthermore, with our cut-off scheme the relative value 
of the anisotropy is rather insensitive to the precise value 

of Cq. 

In Fig. 13 we show the obtained effective interaction 
as a function of the hole concentration -p for two differ- 
ent numbers of fc-points (iV) below the cut-off energy 
E^^tos = 3273K for a Mn separation of i? = a/2 « 2.8A. 
For these calculations we used a dimensionlcss coupling, 
ipd = ^^H^r = 0.2, with S = 5/2 the Mn spin, m the 
free electron mass, and a — 5.65A the lattice constant. 
This value of jpd roughly corresponds to the exchange 
coupling used in Ref. [231 ■ The accuracy of the calcu- 
lation can be inferred both from the difference between 
the N = 3220 and N = 1688 data, and from the size of 
the discrete jumps. It is around 10 — 20% at best. Thus, 
one cannot give a quantitatively reliable estimation of 
the anisotropy from these numerical calculations. In this 
regard it is somewhat misleading that for a fixed num- 
ber of holes the exchange energy behaves very nicely as 
a function of the angle of the two magnetic impuritieS)2S 
as demonstrated in Fig. O In fact, our experience with 
other local density of states calculations shows that one 
may need around 10^ k-points to achieve an accuracy of 
~ 1%, which seems to be beyond the scope of present-day 
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numerics^* 

It is therefore questionable how well these data can 
be trusted in regard to the estimation of the anisotropy 
energy. However, a qualitative assessment can be made. 
While anisotropy energies 8 J are definitely smaller than 
the numerical errors of the total exchange energies, they 
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FIG. 2: Effective exchange J and anisotropy 5, J for two Mn 
ions. We used the cut-off scheme H2()|l with ao = A A, and 
a cut-off energy iScutoff = 3273 K. The main parts show the 
obtained interaction for L = 23a and L — 29a at concentra- 
tions corresponding to "closed sheUs", while the insets show 
the L = 29a results for all concentrations. 



FIG. 3: Spin orientation dependence of the ground state 
energy for a fixed number of holes corresponding to a hole 
concentration of p = 0.2nm~^: (a) One of the spins is aligned 
along the z-axis, and the other is rotated in the xz plane, (b) 
The spins are parallel and rotated in the xz plane. 



still show a consistent pattern (see Fig. |2Jl, suggesting 
that the order of magnitude of the anisotropy SJ/J\\ is 
probably correctly given, even for these small values of 
N. 

The anisotropy in Fig. [21 for a Mn separation of i? « 
2.8A is much smaller than the one obtained within the 
spherical approximation. This somewhat surprising re- 
sult can be understood as follows: for short Mn-Mn sep- 
arations the exchange interaction is dominated by high 
energy electron-hole excitations. However, as shown in 
Fig. EI the four band model (Aso —^ oo) provides a good 
approximation to the exact eigenstates only up to en- 
ergies E ^ 0.15 — 0.2eV, corresponding to hole concen- 
trations up to p '^ 0.3 — 0.4nm^'^. For these high-energy 
excitations it therefore fails and largely overestimates the 
strength of spin-orbit interaction. 

On the other hand, for hole concentrations less than 
p '^ 0.3 nm~'^ the Fermi energy is in a range where 
the A50 — > 00 approximation for the single particle 
states is appropriate. While for short Mn-Mn separa- 
tions, electron- hole excitations at all energy scales con- 
tribute to the RKKY interaction, for larger values of R 
the RKKY interaction is dominated by electron-hole ex- 
citations in the vicinity Sk ^ 1/R of the Fermi surface. 
We expect therefore that the anisotropy will increase for 
these concentrations with increasing Mn-Mn separation 
R. 



These expectations are indeed supported by our nu- 
merical results shown in Fig. [Sl where we find that the 
typical anisotropy, SJ/J\\ increases with the Mn-Mn sep- 
aration and becomes of the order of ~ 20% for the typical 
Mn-Mn separation R ^ lOA at a; = 0.05 Mn content. 

In Fig. El we show the exchange interaction and the 
anisotropy as a function of the separation of the two Mn 
ions for various hole concentrations. The exchange inter- 
action increases monotonically as the hole concentration 
is increased, while the "interaction range" decreases due 
to the increase of the Fermi momentum. The maximum 
of the anisotropy term, on the other hand, decreases for 
larger hole concentrations, since then the Fermi energy 
moves into a range where spin-orbit coupling is of lesser 
importance. Depending on the specific value of hole con- 
centration, the size of the anisotropy is in the range of 
10 — 30% for typical Mn separations. 

Regarding the directional anisotropy (from the cubic 
symmetry of the lattice) mentioned in Ref . 39] we remark 
that while spin anisotropy can induce a frustrated ground 
state and may thus also change the universality class of 
the ferromagnetic transition, directional anisotropy only 
increases the already existing disorder somewhat, and 
does not induce any frustration. As a consequence, it 
does not change the magnetic properties of the system, 
and plays an unimportant role in an already disordered 
ferromagnet. 
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FIG. 4: Comparison of the hole dispersion within the full 
six band model, the four band model, and the spherical ap- 
proximation. The latter two are reasonable approximations 
for Ef < 0.15 — 0.2eV, which corresponds to a hole density 
ofp<0.25-0.4nm-^ 
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FIG. 5: Average anisotropy for 0.2nm~'' < p < 0.3nm~'^. 
We used the same parameters as in Fig. |21 The anisotropy 
increases with distance and becomes of the size ~ 20% for 
R = 2a, which is around the typical Mn-Mn distance for 
X = 5%, indicated by the arrow. 



In the calculations above we neglected disordered 
Coulomb scattering of the valence holes on static impu- 
rities (As antisites and Mn core charges etc.). This type 
of disorder destroys the coherent propagation of the elec- 
trons and therefore reduces the value of both Jn and 8, J 
for separations larger than the electronic mean free path 
i. For the metallic samples i is of the order of the typi- 
cal Mn-Mn separation, and therefore we expect that the 
suppression is not dramatic. While it is not clear how 
much the ratio 5J/ J\\ is reduced by static disorder, it is 
likely that random anisotropy effects are important for 
the disordered samples too. 
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FIG. 6: Position- and hole concentration-dependence of the 
exchange and anisotropy of the RKKY interaction within the 
six band model. 



IV. MONTE CARLO STUDY 

Having obtained the effective kernels, in this section 
we report the results of classical Monte Carlo (MC) cal- 
culations in the spin degrees of freedom used to study the 
implications of anisotropy on the magnetic properties of 
Gai_a;Mn2:As. Similar Monte Carlo studies of GaMnAs 
have also been carried out on other models^^iSLSi. 

The properties of the obtained ferromagnetic state are 
expected to depend on correlations between the positions 
of Mn spins. As we mentioned in the introduction, it is 
presently unknown from the available experimental data 
to what extent the substitutional, magnetically active, 
Mn positions are correlated. It is therefore worthwhile 
to investigate theoretically the trends in magnetic prop- 
erties as a function of increasing correlations in the Mn 
ion positions. We applied the spin-spin interaction ker- 
nels computed in the first part of this paper to moni- 
tor how the magnetic properties such as the saturation 
magnetization, the Curie temperature, the shape of the 
magnetization curves and the spin distributions change 
as Mn positional correlations are gradually built in via 
the procedure described below. 

To simulate positional correlations and to control dis- 
order we first generate a completely random Mn distri- 
bution within the Ga sublattice of an L x L x L cube 
of Gai-ajMua^As, and then relax the Mn ions through 
a standard T = Monte Carlo procedure with nearest 
and next nearest neighbor hopping, assuming a screened 
Coulomb interaction between the Mn ions. It is ex- 
tremely important to use periodic boundary conditions 
in the course of this relaxation procedure, otherwise the 



Mn ions accumulate on the surface of the cube. 

The typical evolution of the average nearest neighbor 
distance with MC time, tMC: is shown in Fig. [3 for two 
different sample sizes and x — 0.05. Correlations are 
formed within a MC time span of tMC ~ 5. A careful 
investigation of these correlations for dilute samples with 
X ~ 0.01 — 0.02 shows that the Mn ions tend to form 
a somewhat distorted BCC lattice with point defects. In 
the following, we shall use this Monte Carlo time as a 
parameter to control the amount of positional disorder. 
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FIG. 7: Average nearest neighbor Mn distance vs. Monte 
Carlo time for L = Vla and L = 17a. 



In all simulations to be described below the Mn 
spins were replaced by their classical angular variables 
S -^S O, which is expected to be a reasonable ap- 
proximation for S = 5/2. To take account of the fi- 
nite mean free path I « 7A of the valence holesSSiS, 
we used an exponential cutoff for the RKKY interac- 
tion: i^para/pcrp(R-) ^ ^para/perp(R-)e"^/'. We alsO 

introduced a hard cutoff for the effective interaction, 
^hard = 2rMn, whcrc TMn IS related to the Mn concen- 
tration, riMn — [47r ?'Mn/3]~^- This hard cut-off has been 
defined to take into account that the RKKY approxima- 
tion does not make any sense beyond the first "shell" of 
neighbors, since Mn ions are very strong potential scat- 
terers. 



A. Monte Carlo results within the spherical 
approximation 

Within the spherical approximation, the results de- 
pend almost exclusively on the hole fraction / and do not 
depend too much on the specific value of the Mn concen- 
tration X. This is because the effective kernel depends 
on the Mn-Mn distance through y — kp^hR- The typical 
values of y depend on /. However, lattice-specific effects 
play a less important role in the dilute limit, where char- 
acteristic distances are typically much larger than the 
lattice constant. 

In Fig.|Slwe show the temperature-dependence of the 
magnetization M = \{fli)\, for different amounts of dis- 



order, as a function of temperature for two different hole 
fractions, / = 0.1 and / = 0.25. A spontaneous mag- 
netization develops at low temperatures in both cases. 
For small Monte Carlo times (large disorder) the tran- 
sition between the paramagnetic and magnetic phases 
takes place rather smoothly, and then the magnetization 
increases approximately linearly with decreasing tem- 
perature, qualitatively similar to many experiments&S 
The Curie temperature (estimated by where the curves 
would intersect the temperature axis if the high tem- 
perature tails are ignored) decreases with decreasing 
disorder, i.e., disorder tends to enhance the transition 
temperaturei22iiii^ While the M{T) curves of the un- 
relaxed samples do not quite look like usual mean-field 
magnetization curves, they become more and more mean- 
field-like upon relaxing the Mn impurity positions. All 
these properties are characteristic of strongly disordered 
magnets and have been reported earlier j^S 

For both hole fractions we find that the magnetization 
tends to a value at T = that is smaller than that of 
a fully polarized ferromagnet. This effect is mostly due 
to anisotropy induced frustration: We recover the fully 
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FIG. 8: Magnetization vs. temperature for hole fractions 
/ = 0.10 (top) and / = 0.25 (bottom) as a function of Monte 
Carlo time. Note that for / = 0.1 the isotropic kernel results 
in nearly full polarization even for the completely disordered 
sample demonstrating the the reduction of the magnetization 
is dominantly due to anisotropy effects. 
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polarized state when we substitute the kernel with its 
angle averaged value. 

Correlations between the Mn sites decrease frustration 
effects and tend to increase this remnant magnetization, 
and for / = 0.1 the fully relaxed system recovers 90% 
of the magnetization (see Fig. I10|l . The corresponding 
evolution of the distribution of the ground state spin ori- 
entations is shown in Fig.O The angle 9 in Fig.j^ldenotes 
the angle with respect to the direction of the ground state 
magnetization vector, n : cos^i = J7i • n. When all spins 
are aligned, P{cos9) = (5(1 — cos6), and the more ordered 
the positions of the Mn ions, the more peaked the distri- 
bution becomes around cos 9^1. 
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FIG. 9: Ground state spin distribution for / = 0.10. The 
angle 9 is measured with respect to the ground state magne- 
tization. 



The simulations with / = 0.25 show a much stronger 
reduction of the zero temperature magnetization relative 
to the case of / = 0.10, even for the relaxed samples. 
Furthermore, for a fixed disorder the Curie temperature 
is reduced by about 20 — 40% for / — 0.25 with respect 
to that of / = 0.1, even if we take into account the factor 
of ~ 3 increase in the energy scale Ej. In view of the 
results of the following section, we believe that these re- 
sults are artifacts of the spherical approximation: As we 
emphasized earlier, within the spherical approximation 
oscillations appear in the interaction kernel, which show 
rather specific features around y = kphR ~ 2, which 
is just the typical value of y for / = 0.25. While for 
/ — 0.10 the parallel component of the kernel is ferro- 
magnetic for typical Mn-Mn separations, for / = 0.25 it 
becomes antiferromagnetic. Therefore it is likely that for 
/ = 0.25 both the anisotropy and the antiferromagnetic 
part of the RKKY coupling play an important role.^^ 

The main effects due to correlations between the Mn 
impurities are summarized in Fig. 1101 where we show 
the T = magnetization for L/a — 17 and the L = oo 
extrapolated value of Tc as a function of Monte Carlo 
time. This latter has been obtained by measuring the 
maximum of the susceptibility for various system sizes 
and then extrapolating to L = oo by using the criti- 
cal exponents for the Heisenberg model known from e 



expansion. ^^ Both quantities change monotonically with 
disorder, with a time scale similar to the one with which 
the disorder changes (see Fig.O. 
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FIG. 10: Top panel: the L — oo extrapolated Curie tempera- 
ture as a function of Monte Carlo time. The dashed curves are 
a guide to the eye; Inset to top panel: The scaling of the peak 
in susceptibility with system size L, for t^fc ~ 2.5;Bottom 
panel: M{T = 0) for / = 0.10 and L = 17o. 



B. Monte Carlo results within the six band model 

In this subsection we study the magnetic properties of 
Gai_a;Mn2;As using the effective interaction kernel com- 
puted within the six band model. In this case the in- 
teraction kernel depends on the specific direction of the 
two Mn ions, and, in principle, one should compute it for 
all possible positions of the two Mn ions, as the effective 
interaction is, in general, a rather complicated tensor. 
This is next to impossible, and besides, we do not expect 
to obtain quantitatively correct results from it anyway, 
since in these calculations we neglect the polarization of 
the holes in the ferromagnetic state, the strong poten- 
tial scattering off the Mn cores, and the effects of various 
defects in Gai_a;Mna;As. 

Instead, to obtain a qualitative picture, we will pur- 
sue the following strategy: By tetragonal symmetry, the 
effective Mn spin-spin interaction kernel reduces to the 
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simple form Eq. (|11() provided that the two Mn ions are 
aligned along the x, y, or z direction. We will there- 
fore use the effective interaction Eq. 1)11(1 . but we will 
substitute the kernels in this expression by the ones we 
computed in Section llTTI In this way we obtain a qualita- 
tively correct description of the interaction between the 
Mn spins in Gai_a;Mn2; As, which captures approximately 
the spin-orbit coupling induced anisotropy effects. 

In this case the structure of the M{T) curves depends 
not only on the hole fraction /, but also on the Mn con- 
centration. In fact, our results show that for the active 
Mn concentration range 0.03 < a; < 0.05, Tc is very sen- 
sitive to the Mn concentration x, but exhibits a much 
weaker dependence on the hole fraction. This originates 
from the specific property of the interaction kernel shown 
in Fig. 1^ As we showed in the previous section, the role 
of spin-orbit coupling induced anisotropy is also more 
pronounced for larger Mn-Mn separations, i.e., smaller 
Mn concentrations. 
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FIG. 11: Temperature-dependence of the magnetization for a 
sample of linear size L = Q a with active Mn concentration 
X = 0.05, hole fractions / = 0.25 and / = 0.4 as a function 
of Monte Carlo relaxation time. The insets show the low 
temperature saturation of the magnetization. 



between the ground state of the sample and an individual 
spin of the order of ^ 16 degrees. This non-coUinearity 
disappears once we introduce correlations between the 
Mn sites. These results clearly show that for larger Mn 
concentrations the spherical approximation badly fails 
and a more complete six band model must be used. 

We also find that for even smaller concentrations, 
X < 0.03, the spin-orbit coupling induced disorder plays 
a more important role, and the obtained ground state 
is non-collinear, similar to the one obtained within the 
spherical approximation. These results indicate that 
non-collinear states may appear for smaller Mn concen- 
trations. This small concentration regime is, however, 
definitely out of reach for an RKKY approach: At these 
small concentration disorder also plays an important role 
and most likely an impurity band description of the ma- 
terial is necessarjsiSiSfliil. 
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FIG. 12: Temperature-dependence of the magnetization for a 
sample of linear size L = 10 a with active Mn concentration 
X — 0.03, hole fractions / — 0.25 and / — 0.4 as a function 
of Monte Carlo relaxation time. The insets show the low 
temperature saturation of the magnetization. 



Typical magnetization curves are shown in Figs. 1111 
and [la for a; = 0.05 and x = 0.03. For x = 0.05 the 
ground state is always almost fully polarized, even for 
the fully disordered, unrelaxed sample, only Tc is gradu- 
ally suppressed upon relaxing the impurity positions. For 
X — 0.03, however, a slight non-coUinearity appears for 
the fully disordered sample, corresponding to a 4% sup- 
pression of the total magnetization and a typical angle 



V. CONCLUSIONS 

In this paper we continued the route of Ref. [2^ to 
build effective spin models for metallic Gai_a;Mna;As. 
We constructed interactions between the Mn spins medi- 
ated via indirect exchange by using two different valence 
band models. First we computed the effective interaction 
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within the spherical approximation, where the calcula- 
tions can be performed analytically, and then we studied 
the effective interaction numerically using the more com- 
plete six band model, also studied in Ref. |39| . 

We find, in agreement with Ref. [S^l , that the spherical 
approximation badly fails for short Mn separations where 
it overestimates the strength of the spin-orbit coupling in- 
duced anisotropy by more than one order of magnitude. 
However, we find using the six band model, that the 
strength of the spin-orbit coupling induced anisotropy in- 
creases with Mn-Mn distance, and becomes of the order 
of ~ 20%, qualitatively (but not quantitatively) agreeing 
with the results obtained using the spherical approxima- 
tion and also in rough agreement with the experimentally 
observed FMR linewidthi^LS This implies that the frus- 
tration effects discussed in Ref. [23 are probably more 
pronounced for small Mn concentrations. 

We carried out classical Monte Carlo simulations to 
study the implications of the computed effective inter- 
actions. Within the spherical approximation, we always 
find a generically non-collinear state due to orientational 
frustration (random anisotropy). It has been speculated 
earlier that this state may be a ferromagnetic state of 
spin glass nature. ^^ However, our studies of hysteresis 
indicate that this random anisotropy results in a con- 
ventional, though non-collinear and disordered ferromag- 
netic state. This result seems to be supported by the fact 
that random anisotropy is presumably irrelevant at the 
Heisenberg fixed point, ^® and therefore does not change 
the critical scaling of a Heisenberg magnet. 

The strongly non-collinear states found within the 
spherical model are partly an artifact of the spherical 
approximation, which is only valid at sufficiently small 
hole concentrations. Using the effective interaction ob- 
tained within the six band model of the valence holes, we 
find that the anisotropy effects are too weak to induce a 
strongly non-collinear state for (active) Mn concentra- 
tions as large as x = 0.05. However, anisotropy becomes 
more important for lower concentrations, where it may 
induce a frustrated and non-collinear stateiiS In partic- 
ular, for X = 0.03 we find that the ground state of the 
fully disordered sample is not fully coUinear, and indi- 
vidual spins deviate from the ferromagnetic orientation 
of the sample by about ~ 16 degrees. This tendency is 
expected to be even more pronounced for lower Mn con- 
centrations where we expect a non-collinear ferromag- 
netic state. ^^ Indeed, for Mn concentrations in the range 
X = 0.015 — 0.02 we find a non-collinear state similar to 
the one obtained in the spherical approximation, though 
this concentration range may already be well out of the 
range of the metallic approximation used in this p aper, 
and the impurity band approach of Refs. |l9l2Cl4ll | may 
be more appropriate in this limit. Experimental evidence 
also supports the presence of a non-collinear state in the 
regime of small Mn concentrations!^ 

Let us emphasize again, that the Mn concentrations 
above denote the concentration of active substitutional 
Mn ions, which can be substantially smaller than the 



nominal Mn concentration of the samples. ^'^ In particu- 
lar, it is known that a substantial amount of the dopant 
Mn ions go to interstitial sites in course of the out-of 
equilibrium growth processjSiS and these interstitial Mn 
ions are also believed to bind to the substitutional Mn 
ions and "neutralize" them from the point of view of 
ferromagnetism.^^ In this way, it is quite possible that 
an unannealed sample with a nominal Mn concentration 
of Xnom '^ 0.05 has an active Mn concentration of only 
around x — 0.01 — 0.03 and have a non-collinear ground 
state. It is, on the other hand quite likely that annealed 
samples and samples with higher Mn concentrations form 
a ferromagnetic state where the Mn spins are almost per- 
fectly coUinear. 

As first pointed out in Ref. [23, the observed satu- 
ration magnetization of Gai_a:Mn2,As is much less even 
in annealed samples than expected based on the nomi- 
nal Mn concentration of the samples. The orientational 
frustration discussed in this paper provides a possible 
mechanism to explain the missing magnetization in unan- 
nealed and "underdoped" samples, where one can sub- 
stantially increase the saturation value of the magne- 
tization, M{T — 0) with a relatively small magnetic 
field compared to Tc, H < 1 T^ Other mechanisms 
have also been proposed: The main mechanism that re- 
duces the magnetization of annealed samples seems to be 
provided by the presence of interstitial Mn ions, which 
- as mentioned earlier - "neutralize" substitutional Mn 
impuritiesii^ It has also been proposed that As antisite 
defects can induce a non-collinear statej^ and an intrin- 
sic non-coUinearity has been proposed also in Ref. (26i27l| . 
Finally, quantum fluctuations due to the antiferromag- 
netic coupling between the Mn spins and the valence 
holes have been also proposed recently. ^^ It is probably a 
combination of all these effects that is finally responsible 
for the missing magnetization of Gai-ajMn^As. 

Finally, we comment on the approximations made in 
this paper. Throughout, we followed a common practice 
in the literature^i and neglected the potential scattering 
off the Mn ions. This approximation is expected to pro- 
duce qualitatively reliable results in the metallic limit con- 
sidered and we believe that all the trends reported here 
are robust. We also neglected the polarization of the va- 
lence holes induced by the ferromagnetically alligned Mn 
spins, which may be important at low temperatures. In 
summary, we believe the answers to the questions of how 
anisotropic the Mn spin-spin interactions are in metallic 
GaMnAs and how anisotropy interplays with Mn posi- 
tional correlations to affect important magnetic proper- 
ties are now on more solid ground, though further studies 
are needed to understand frustration effects in the regime 
X < 0.03. 

NOTE: While preparing this work for publication, a 
preprint appeared^ where the authors arrive at rather 
similar conclusions using a very different, tight binding 
approach. 
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APPENDIX A: 



CALCULATION OF RKKY 
KERNEL 



In this appendix, we derive explicit expressions for the 
kernels Xpar(i?) and Kpcip{R) appearing in Eq. ©. By 
spherical symmetry, we can assume that the two Mn im- 
purities are alligned along the z axis. In this case, the 
kernel K""^ becomes diagonal, K"^ = K"Sai3, and can 
be written by converting the sums in Eq. (UJ to integrals 
and assuming a parabolic dispersion for the light and 
heavy hole bands as 
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27r2 



k'^dk' 




1 




27r2 


fc'2 

2m„, 


- 


fc2 
2m„ 



(I J" , (e, e', 0, ,^')l' 2 cos(fci?cos(0) - k'Rcos{0')))e,e',4>,4>' 



r 



iM 



Here we have taken the angles 9 and cj) 
(0' and (j)') parametrize the directions k = 
(sin 6*008 0, sin 6* sin 0, cosO) and k', and the brackets 
denote the angular average over the angles 0,9',(f) and 
(/)'. Making the substitutions k — > ^J2m^€Fk and 
k' —>■ ^/2rnfjepq yields the expression 



if"(R) = 8eFVg^V I dk ( dq 
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where 
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and we introduced the dimensionless couplings, g^ — 

711 1/2 

To compute K^^CR) we first evaluate 
I^^,{kp,jj_kR,kF,fj.'qR) by performing the angular 
integrals. Since lan'i^F.fikR^kF^fj^iqR) only depends on 
k and q through the combinations kp^^kR and kp^^iqR, 
it is worth defining the heavy hole contributions to it as: 
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The contributions /^;, /j'^, and Ifl, and the corresponding 
contributions to the kernel, K'j^^^, K"f^, K"i, K^, can be 
defined in an analogous way. 

We demonstrate the procedure of computing K" by 
the example of the heavy hole contribution to K^^f^. It is 



straightforward to evaluate the angular integrals, and for 
I^f^ one obtains the following expression: 
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where F and G have the obvious definitions. The heavy 
hole contribution to K^{R) thus reads 
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As shown in Appendix IbI for F and G of the form that 
appear in Eq. I|A2|1 the integrals over q in Eq. I|A7|1 can 
be carried out to give 



qk 



*A '%^~k^ 
1 



F{k)F{q) 



TT 



dkk 



F{k)Ymx 



q-tO 



F+(kf -F~{kf 
2g2 



k^ 



F+{q) 



(A8) 



14 



where F^ and F~ denote the parts of the function 
F = F+ + F^ which are analytical on the upper and 
lower half-planes, respectively. Applying this formula, 
we obtain upon making the substitution a = kphR and 
h = kp^iR 
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It is also possible to evalute the integral IJA9I) . but the 
resulting formula is very long, and for practical purposes 
it is simpler to evaluate Eq. l|A9p numerically. 

The remaining parts of the kernel can be evaluated 
in a similar way, and are given below. By symmetry, 
K^ = Ry for this arrangement, and therefore only the 
x-component is displayed. The spatial dependence of the 
various parts of the kernel is shown in Fig. ^ 
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APPENDIX B: EVALUATION OF SINGULAR 
RKKY INTEGRALS FOR SPIN 3/2 PARTICLES 



In this appendix we establish the identity Eq. 
First we note that the functions F{k) that appear in 
the evaluation of the RKKY kernel, Eq. ljA2p . have two 
important properties that we will use in course of the 
derivation: (i) they are odd, F{—k) = —F{k) and (ii) F 
is regular at the origin. 

First let us prove that 



lim5_ 



1+1(5 



dz 



l+iS 



dq- 



qz 



q- 



;F{z)F{q) = 0, 



(Bl) 
where 5 is a positive infinitesimal quantity. Decomposing 

the integrand i'' .> = ^ ( — — I — ^ ) , introducing the 
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variable x — z ~ i6 I, and using the identity 
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where P denotes the principal part of the integral. The 
principal value integral over the first term is zero by sym- 
metry, and the integral over the last two terms terms sim- 
ply evaluates to i^{qF{q)F{q) - {-~q)F{~q)F{q)) = 0, 
since F{—q) = —F{q). 

We now return to the formula, Eq. (|A8|) . Using 
Eq. HBll) and the fact that the integrand in Eq. (|A8|) 
is even in both k and q, we can extend the region of 
integration to obtain 
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FIG. 13: The contours used in the evaluation of Eq. I|B4 
The dots indicate possible singularities of the integrand. 



where C denotes the contour z E [—l + id,l + iS]. 
We now turn our attention to the integral 



R 



dq 



1 



1 



q — z q + z 
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We first deform the contour of q such that it passes below 
the origin using the fact that the integrand of Eq. IjBSp 
is analytical at the origin. Then we decompose the func- 
tion F as i^ = F+ + F~ , with F+ and F~ being analyt- 
ical functions in the upper and lower half-planes (apart 
from possible singularities at the origin). For example, in 
Eq. lli, we rewrite F{q) as Fiq) = sM ^ g-^^^, 
and we can similarly decompose G{q) in Eq. lASp in the 
same way term by term. We can now close the contours 
for the part ~ F~^ on the upper half-plane, and similarly 
for ^ F~ on the lower half-plane as shown in Fig. E| to 
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obtain 
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with z ~ k + iS. Evaluation of the first integral 



and the sec- 
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ond integral gives 27ri(— 1)F^(— fc). Using the property 
F^{—k) ~ — F+(fc) we obtain the final result 
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Since the worst singularity is F^ (q) ^ q ^ , we can write 
the residue part of this expression as 
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and we finally obtain 

j-ldkj-^dq^F{k)F{q)= (B7) 

i j-\ dkkF(k)2m [2F+(fc) + hniq^o {^F+{q)^' 



In the first term we can make the replacement under 
the integral F+{k) = i (F+(fc) + F+(-fc)] since kF{k) 
is even. Recalling that F{q) ~ F^{q) + F^{q) and 
F+(— fc) = —F^{k) we then obtain, 



''h'^^ 

1 



F{k)F{q) 



(B8) 



■n 



= - / dkk 







F+{kf -F-{kf 



-^(fc)limq^O , ^2_fc2^ 

which is just the identity we wanted to establish. 



(B9) 



^ For recent reviews see J. Konig et al. in Electronic Struc- 
ture and Magnetism of Complex Materials, edited by 
D.J. Singh and D.A. Papaconstantopoulos (Springer Ver- 
lag 2002); R. N. Bhatt et al, J. Superconductivity INM 
15, 71 (2002); F. Matsukura and H. Ohno and T. Dietl 
in:Handbook on Magnetic Materials, Elsevier 2002. 

2 D. P. DiVincenzo, Science 270, 255 (1995). 

^ M. A. Nielson and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, UK, 2000). 

* D. D. Awschalom, N. Samarth, and D. Loss, eds.. Semicon- 
ductor Spintronics and Quantum Computation (Springer 
Verlag, Heidelberg, 2002). 

^ I. Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 76, 
323 (2004). 

^ K. W. Edmonds, P. Boguslawski, K. Y. Wang, R. P. Cam- 
pion, N. R. S. Farley, B. Gallagher, C. T. Foxon, M. Saw- 
icki, T. Dietl, M. B. Nardelli, et al., Phys. Rev. Lett. 92, 
037201 (2004). 

'' H. Ohno, S. Shen, F. Matsukura, A. Oiwa, A. Endo, 
S. Katsumoto, and Y. lye, Appl. Phys. Lett. 69, 363 
(1996). 

* S. J. Potashnik, K. C. Ku, S. H. Chun, J. J. Berry, 
N. Samarth, and P. Schiffer, Appl. Phys. Lett. 79, 1495 
(2001). 

® K. M. Yu, W. Walukiewicz, T. Wojtowicz, I. Kuryliszn, 

X. Liu, S. Sasaki, and J. K. Furdyna, Phys. Rev. B 65, 

201303 (2002). 
^'^ T. Hayashi, Y. Hashimoto, S. Katsumoto, and Y. lye, 

Appl. Phys. Lett. 78, 1691 (2001). 
" K. M. Yu, W. Walukiewicz, T. Wojtowicz, W. L. Lim, 

X. Liu, Y. Sasaki, M. Dobrowolska, and J. K. Furdyna, 

Appl. Phys. Lett. 81, 844 (2002). 
^^ B. S. Sorenson, P. E. Lindelof, J. Sadowski, R. Mathieu, 

and P. Svedlindh, Appl. Phys. Lett. 82, 2287 (2003). 
^^ H. Ohno, J. Magn. Magn. Mater. 200, 110 (1999). 



" K. W. Edmonds, K. Y. Wang, R. P. Campion, A. C. Neu- 
mann, N. R. S. Farley, B. L. Gallagher, and C. T. Foxon, 
Appl. Phys. Lett. 81, 4991 (2002). 

15 K. C. Ku, S. J. Potashnik, R. F. Wang, S. H. Chun, 
P. Schiffer, N. Samarth, M. J. Seong, A. Mascarenhas, 
E. Johnston-Halperin, R. C. Meyers, et al., Appl. Phys. 
Lett. 82, 2302 (2003). 

^^ S. J. Potashnik, K. C. Ku, R. F. Wang, M. B. Stone, 
N. Samarth, P. Schiffer, and S. H. Chun, J. Appl. Phys. 
93, 6784 (2003). 

^■^ S. J. Potashnik, K. C. Ku, , R. Mahendiran, S. H. Chun, 
R. F. Wang, N. Samarth, and P. Schiffer, Phys. Rev. B 66, 
021408 (2002). 

1* L. Bergqvist, P. A. Korzhavyi, B. Sanyal, S. Mirbt, L A. 
Abrikosov, L. Nordstrom, E. A. Smirnova, P. Mohn, 
P. Svedlindh, and O. Eriksson, Phys. Rev. B 67, 205201 
(2003). 

1^ G. A. Fiete, G. Zarand, and K. Damle, Phys. Rev. Lett. 
91, 097202 (2003). 

2° M. Berciu and R. N. Bhatt, Phys. Rev. Lett. 87, 107203 
(2001). 

•^1 For a recent topical review on defects in GaMnAs see: C. 
Timm, J. Phys. Cond. Matt. 15, R1865 (2003). 

^2 J. Bhnowski and P. Kacman, Phys. Rev. B 67, 121204(R) 
(2003). 

^^ P. A. Korzhavyi, L A. Abrikosov, E. A. Smirnova, 
L. Bergqvist, P. Mohn, R. Mathieu, P. Svedlindh, J. Sad- 
owski, E. L Isaev, and Y. K. V. O. Eriksson, Phys. Rev. 
Lett. 88, 187202 (2002). 

2" F. Maca and J. Masek, Phys. Rev. B 65, 235209 (2002). 

2''' G. Zarand and B. Janko, Phys. Rev. Lett. 89, 047201 
(2002). 

'^^ J. Schliemann, Phys. Rev. B 67, 045202 (2003). 

^'^ J. Schliemann and A. H. MacDonald, Phys. Rev. Lett. 88, 
137201 (2002). 

2* J. D. J. Priour, E. H. Hwang, and S. D. Sarma, Phys. Rev. 



17 



Lett. 92, 117201 (2004). 
^^ A. Singh, A. Datta, S. K. Das, and V. A. Singh, Phys. Rev. 

B 68, 235208 (2003). 
^° G. Bouzerar, J. Kudrnovsky, and P. Bruno, Phys. Rev. B 

68, 205311 (2003). 
^^ C. Timm and A. H. MacDonald, cond-mat/0405484 
^^ J. Konig, H.-H. Lin, and A. H. MacDonald, Phys. Rev. 

Lett. 84, 5628 (2000). 
^^ M. Abolfath, T. Jungwirth, J. Brum, and A. H. MacDon- 
ald, Phys. Rev. B 63, 054418 (2001). 
^'' C. Timm, F. Schafer, and F. von Oppen, Phys. Rev. Lett. 

89, 137201 (2002). 
^^ P. Redlinski, (private communication). 
^^ A. Baldereschi and N. O. Lipari, Phys. Rev. B 62, 2697 

(1973). 
^'^ J. Konig, T. Jungwirth, and A. H. MacDonald, Phys. Rev. 

B 64, 184423 (2001). 
^* W. Kohn and J. M. Luttinger, Phys. Rev. 98, 915 (1955). 
^^ L. Brey and G. Gomez-Santos, Phys. Rev. B 68, 115206 

(2003). 
*° A. C. Hewson, The Kondo Problem to Heavy Fermions 

(Cambridge University Press, Cambridge, England, 1997). 
"^ M. Berciu and R. N. Bhatt, Phys. Rev. B 69, 045202 

(2004). 
"^ J. S. Blakemore, J. Appl. Phys. 53, R123 (1982). 
*^ M. Linnarson, E. Janzen, B. Monemar, M. Kleverman, and 

A. Thilderkvist, Phys. Rev. B 55, 6938 (1997). 
"** J. Szczytko, A. Twardowski, K. Swiatek, M. Palczewska, 

M. Tanaka, T. Hayashi, and K. Ando, Phys. Rev. B 60, 

8304 (1999). 
*^ J. Okabayashi, A. Kimura, O. Rader, T. Mizokawa, A. Fu- 



jimori, T. Hayashi, and M. Tanaka, Phys. Rev. B 58, 
R4211 (1998). 

*^ K. Yosida, Theory of Magnetism ( Springer- Verlag, Berlin, 
1996). 

*^ T. G. Rappoport, P. Redlinski, X. Liu, G. Zarand, J. Fur- 
dyna, and B. Janko, Phys. Rev. B 69, 125213 (2004). 

** A similar, but more complicated expression gives the struc- 
ture of Jpd(k, k') in GaMnAs, too. 

*^ The time required for the numerical calculation scales as 

^'^ J. Schliemann, J. Konig, and A. H. MacDonald, Phys. Rev. 

B 64, 165201 (2001). 
^^ M. P. Kennett, M. Berciu, and R. N. Bhatt, Phys. Rev. B 

66, 045207 (2002). 
^^ H. Ohno, H. Munekata, S. von Molnar, and L. L. Chang, 

Phys. Rev. Lett. 68, 2664 (1992). 
^^ A. V. Esch, L. V. Bockstal, J. D. Boeck, G. Verbanck, A. S. 

van Steenbergen, and P. J. Wellmann, Phys. Rev. B 56, 

13103 (1997). 
^'^ R. N. Bhatt and X. Wan, Int. J. Mod. Phys. C 10, 1459 

(1999). 
'^^ C. Zhou, M. P. Kennett, X. Wan, M. Berciu, and R. N. 

Bhatt, Phys. Rev. B 269, 144419 (2004). 
^^ J. C. L. Guillou and J. Zinn- Justin, Phys. Rev. B 21, 3976 

(1980). 
" X. Liu, Y. Sasaki, and J. Furdyna, Phys. Rev. B 67, 205204 

(2003). 
'^^ J. T. Chayes, L. Chayes, D. S. Fisher, and T. Spencer, 

Phys. Rev. Lett. 57, 2999 (1986). 
"'^ A. H. MacDonald, (private communication). 



